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Abstract 

Effective field theory (EFT) has been recently used for the calculation of neutron-deuteron radiative capture at very 
low energies. We present here the use of EFT to calculate the two-body photodisintegration of the triton, considering 
the three-body force. The calculated cross section shows sharply rising from threshold to maximum about 0.88 mb at 
~ 13 MeV and decreasing slightly to about 0.81 mb at ~ 19 MeV, in corresponding to the experimental results. Our 
results are in good agreement with the experimental data and the other modern realistic two- and three-nucleon forces 
AV18/UrbanaIX potential calculation. 
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1. Introduction 

The total photodisintegration cross sections of the three-body nuclei are important quantities for understanding 
electromagnetic processes in nuclei. The three-nucleon systems involving nucleon-deuteron radiative capture and 
photodisintegration of 3 He and 3 H have been investigated in experimental and theoretical works over the past decades. 

In the experimental investigations, a few measurements of triton photodisintegration have been done. More re- 
cently, including of the two- and three-body reactions, photodisintegration cross sections are measured at low ener- 
gies 111]. Mono energetic photons are used and neutrons are detected, in this experiment. The obtained 3 H(y,n) 2 H 
cross section shows sharply rising from threshold to the maximum of 0.9 mb at ~ 12 MeV and decreasing only slightly 
to 0.8 mb at ~ 19 MeV. The results are the two-body breakup cross sections for 3 H (]]]. The same shape for the 3 He 
cross section has been reported but with the lower magnitude. 

During the past several years, different theoretical approaches have been used in order to investigate the cross 
sections of the nucleon-deuteron radiative capture and photodisintegration of 3 He and 3 H from the very low energy 
up to 100 MeV. The main aim of these theoretical studies has been the calculations of the bound state and scattering 
wave functions in order to be used in matrix elements of the electromagnetic current operators and ultimately testing 
the underlying two, three body interactions as well as currents with certain degree of complexity for the three-nucleon 
systems. In many theoretical treatments, by solving the three-nucleon Faddeev equations based on simple finite rank 
forces, it is found that the interacting three-nucleon continuum is crucial for understanding the process in three- 
body reactions J2]. Sandhas and his collaborators investigated many theoretical treatments of photo-effect in 3 H 
and 3 He three-nucleon systems by employing the realistic nucleon-nucleon forces exclusive process at 0-40 MeV 
energies 01 ■ For covering a larger energy range, the inclusive two- and three-body breakup cross sections have been 
studied by Efros et al. @]. Recently, the Bonn group used the modern nucleon-nucleon forces. In these works, the 
current was restricted to the dominant E\ multipole Ja] or to the E\ and £2 multipoles 0]. They compared the Faddeev 
approach and a hyperspherical harmonic expansion together with the Lorentz integral transform method using AVI 8 
together with Urbana IX. 



"Corresponding author 

Email address: H-Sadeghi@araku . ac . ir (H. Sadeghi ) 
'Tel:+988614173400, Fax:+988614173406 



Preprint submitted to Few-Body Systems 



April 12, 2010 



For the p-d and n-d radiative capture, a magnetic dipole Ml transition is a dominant contribution at the very low 
energy, in the plane wave(Born)approximation |8|]. In such investigations, the authors employed Faddeev calculations, 
with the inclusion of three-body forces and meson exchange current(MEC). More recently, a very good agreement 
reached with the calculation of the total three-nucleon photodisintegration cross section with two quite different ap- 
proaches |@| . Rather detailed investigation of such processes has been performed by Viviani et al. fioll . Recent 
calculation using manifestly gauge invariant currents reduced the spread Hill , but the result including three-body 
currents, 0.558 mb, still over-predicts the cross-section by 10%. The Vucom potential leads to a very similar descrip- 
tion of the cross section as the AVI 8 interaction, considering the Urbana IX three-body force for photon energies 
45 < a> < 120 MeV, while larger differences are found close to threshold 111 211 . 

Du ring the last few years, nuclear pionless EFT has been applied to the two-, three-nucleon systems 



16l I17L I18L I19L 12111 and recently developed pionless EFT is particularly suited to the high order precision calculation. 



A recent model independent neutron-proton radiative capture calculation has been done, using EFT to predict a theo- 
retical uncertainty 1 % for center of mass energies E < 1 MeV I22I1 . At these energies, which are relevant for big-bang 
nucleosynthesis, isovector magnetic transitions M\ and isovector electric transitions E\ give the dominant contribu- 
tions. In past years, we have calculated the cross section of nd — > 3 fly process II231 l24l 12511 . No new three-nucleon 
forces are needed in order to achieve cutoff independent results of neutron-deuteron radiative capture process up to 
next-to-next-to leading order(N 2 LO). The three body parameters will be fixed by the triton binding energy and Nd 
scattering length in the triton channel. 

In this paper, we study the two-body photodisintegration of the triton y 3 H — > nd using pionless EFT and insertion 
of the three-body force, up to N 2 LO. The evaluated cross section has been compared with experimental data and the 
total cross section of the three-nucleon photodisintegration calculation of the modern realistic two- and three-nucleon 
forces AV 1 8/UrbanaIX potential models. Close agreement between the available experimental data and the calculated 
cross section is reached. We demonstrate convergence and cutoff independence order by order in the low energy 
expansion. 

The paper is organized as follows: in Section [2] we briefly review theoretical framework for calculating the 
Faddeev integral equation, Lagrangians for the reaction involved in our consideration, the relevant diagrams, three- 
body forces and cutoffs dependence. Comparisons of our results with the corresponding experimental and theoretical 
data are given in Section [3] Summary and conclusions follow in Section |4] 



2. Theoretical framework 

In the present work we confine ourselves with the electric and magnetic dipole approximation for the transition 
operators. These transitions are known to make most of the total cross section at least for the low energy region 
where the pionless EFT theory is most reliable. The photodisintegration reaction, y 3 H — > nd, in these energies occur 
predominantly via electric and magnetic, specially electric-dipole transition. The triton photodisintegration cross 
section is then given by fl: 

cr 7 ^nd(E 7 ) = (2n)\—)M(E y )E y , (2. 1) 

where E y is the incident photon energy and M{E y ) will be calculated from electric- and magnetic-dipole transition 
amplitudes. These transition amplitudes have been prepared by considering of the Faddeev equation for final nd- 
scattering and for the initial triton bound state system to some order (e.g. LO, NLO and N 2 LO). We then take these 
Faddeev amplitudes and folded with the photon-interaction with nucleons when the photon kernel is expanded to the 
same order. 

The Lagrangians relevant to photodisintegration of the triton can be divided into: 

(a) Lagrangian that contributes to nucleon-deuteron scattering. The derivation of the integral equation describing 
neutron-deuteron systems has been discussed before, see e.g. ifm 21, 24]. Only the results will be used here. 



Two mixed configuration in the three-nucleon system are exist. The spin-triplet auxiliary deuteron field D, combines, 
depending on their spin directions, with the nucleon N to total quartet channel spin (S = |) or doublet channel^ = I). 
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Figure 1: The Faddeev equation for Wrf-scattering and diagrams for adding photon-interaction to the Faddeev equation contribution to E\ am- 
plitude. Wavy line shows photon and thick line is propagator of the two intermediate auxiliary fields D s and D,; %: propagator of the exchanged 
nucleon, see eq.(2.2). 



The Faddeev equation is introduced as a coupled equation and is given by 12 ill : 
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where E := |S - ^ - ie is the total non-relativistic energy; the amplitude f9 for the A^, — > A^-process; f9 for the 

Nd t — * AWj-process; k the relative momentum of the incoming deuteron; p the off-shell momentum of the outgoing 
one, with p — k the on-shell point; and the deuteron propagator only at LO is: 
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The projected propagator of the exchanged particle on angular momentum I is given by: 



(2.3) 
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where 2;(z), /th Legendre polynomial is: 



<2i(z) 



(0 



(b) Lagrangian that includes the new term generated by the two-derivative three-body force. Contrary to the terms 
without derivatives, there are different, inequivalent three-body force terms with two derivatives. But only one of 
them, H2 is enhanced over its naive dimensional estimate, mandating its inclusion at N 2 LO 0~8, 21]. 

In that case, a unique solution exists in the 2 S 1/2-channel for each A and vanishing three-body force, but no unique 
limit as A — > 00. As long distance, phenomena must however be insensitive to details of the short-distance physics 
(and in particular of the regulator chosen), Bedaque et al. 1 15, 161, 18 , 2(J showed that the system must be stabilized 
by introducing of a three-body force in a form: 



2 °° / 

n=0 * 



ME + yj 
A 2 



(2.5) 



which absorbs all dependence on the cutoff as A — > 00 |l8ll . At leading order, y «; A, the essential observation are 
of order 0(Q/A) and are independent of energy (or momentum) and hence can be made to vanish by only Hq and not 
any of the higher derivative three-body forces. At NLO, there is only a perturbative change from the LO asymptotic 
of order 0(1/ A 2 )(see JHh). 
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For N 2 LO order, a new three-body force seems to be required. The terms of order (Q/A) 3 is seems to be required, 
which guarantees that the amplitude is cutoff independent up to (Q/A) 2 . The terms will be proportional to ME arising 
from expanding the kernel in powers of q/A and y/ME/A and a three-body force term which contains a dependence 
on the external momenta k 2 and ME can absorb them. It means we need H2 for calculations up to N 2 LO. It has been 
shown that the three body system at N 2 LO can be renormalized without the need for an energy dependent three-body 
force at this order luOll . However, this work only discusses three-nucleon scattering at cutoffs, up to 500 MeV, that are 
large compared to the pion mass. 

(c) Lagrangian that describing currents, contributing to nd — > 3 Hy, which are not included in the previous two 
Lagrangians. 

The polarization can be characterized by the three-vector e, which satisfies the Lorentz condition: e ■ k — 0, where 
k is the unit vector along the photon three-momentum. The sum over the photon polarizations is given by: 

2«£M = Sab -£ a £ b , 
1=1,2 

where the upper index i numerates the two possible independent polarization vectors of the photon. 

For photo production processes or radiative capture, the photon can be characterized by the smallest value of j. In 

this case, we can easily write the corresponding combinations of polarization e and unit vector k for the photon states 
with low multipolarity: 

t — * El (electric dipole), 

e x k — > Ml (magnetic dipole), 
e a kb + ebk a = E a h — > E2 (electric quadrupole), 
x k) a kb + (e x k)bk a = M a b — > M2 (magnetic quadrupole). 



These are the basic formulas for the construction of the matrix elements for electromagnetic processes ( in near 
threshold conditions). Main part of its contribution will be E\ transitions for photon energies below 20 MeV. If the 
ground state of 3 H is thought to be an almost pure 2 S 1/2 state, an electric-dipole transition gives a 2 Pi/2,3/2 final state. 

Cross sections for transitions to final states with isospin T — 3/2 correspond to the three-body breakup while those 
to T — 1 /2 are known to correspond predominantly to the two-body breakup. The resulting matrix element for E\ 
transition of nd — > 3 Hy process, can be written as: 

(*W -&N), (2.6) 

where N, t, q and D are the 2-component spinors of initial nucleon, final 3 H, the unit vector along the 3-momentum 
of the colliding nucleons and the polarization vector of the deuteron, respectively. The spin structure of the matrix 
elements for N + d radiative capture is complicate also for the low energy interaction, as we have here 3 indepen- 

1 + 3 + 

dent multipole transitions: / = - — > Ml and j p — — — > M\ and E%, The following parametrization of the 
corresponding contributions to the matrix element are: 

i(t^N)0 ■ e* x k), 

(fV fl A)(D x [<? x k]) a , (2.7) 

t\&- e* d-£+a--k6-e)N, 

The first two structures in ( 12.7b are correspond to the M\ radiation. 

The E\ and M\ amplitudes receive contributions from the electric and magnetic moments of the nucleon and from 
four-nucleon operators coupling to the magnetic field, which are described by the two-body currents contributing to 



nd -> 3 Hy or 3 Hy -> nd J22J]: 
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£™ = ^T N ^ k ° + kiT 3 )cr.BN + e L B (N T PjN)\N t P 3 N)Bj 

+ l -e L E (N T (T)j P ( V - T)j)N)HN T P a N)Ej + H.C. (2.8) 

where E, B and P a are the electric, magnetic fields and the projector for the 'So channel, respectively. The 
covariant derivative is: 

1 + T-i 

L)K = d lI + ie -A" , 

2 

and P' 1 * are the spin-isospin projectors for the isotriplet, spintriplet channel 

Jf > = -t-^CTi T 2 T 3 , Tr^Pf ] = Itfy. (2.9) 

V8 1 2 

The and are determined from the two-body neutron-proton radiative capture and photodisintegration of the 
deuteron data, respectively, see 1E2I1 for more details. 

The diagrams of E\ contribution is shown in fig. 1. This figure shows simple inclusion of photon to the Faddeev 
equation contribution to £1 amplitude. One of the other possible diagrams which could be contributed in the calcu- 
lation is the insertion photon to exchanged nucleon(middle diagram in fig. 1). This diagram gives a tiny contribution 
when the photon couples electrically to exchanged nucleon, i.e. via P ■ A. This is due to the fact that the initial and 
final state in this process are nearly orthogonal. 

A higher dimension operator coupling an E\ photon to a 3 5i-dibaryon and a p-wave two nucleon final state is 
suppressed by Q 3 in the power counting. This calculation involve only E\ contribution. Meson exchange currents is 
another diagram to be incorporated and it is seen that electric quadrupole E2 contributions are needed to achieve the 
better agreement between theory and experiment, specially at high energies. For higher energies, the present results 
are not complete and calculations with contribution must be carried out. 

The Mi amplitude receives contributions from the magnetic moments of the nucleon and dibaryon operators 
coupling to the magnetic field. Fig. 2 shows the diagrams for adding photon to the Faddeev equation up to N 2 LO. 
There is no interaction of Ho with a photon, because Ho has no derivatives, so it is not affected by the minimal 
substitution P — > P -eA. Contribution from the photon couples to the three-nucleon field (H2) is not negligible(fig. 2) 
and should be calculated up to order of calculation. This coupling comes from the kinetic energy insertion of the 
three-nucleon field. 




Figure 2: The Faddeev equation for Ay-scattering and some diagrams for adding photon-interaction to the Faddeev equation contribution to M\ 
amplitude. The small circles show magnetic photon interaction. Vertices for photon interaction to two- and three-body system have been shown by 
L\ and Hi, respectivly. Remaining notation as in fig-H 
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3. Results and discussion 



The numerical calculations are done using B = 2.225 MeV (yj = 45.7066 MeV)for the AW triplet channel 
(deuteron binding momentum), a s = -23.7 14fm for the AW singlet channel an So scattering length. At LO and NLO, 
Ho is the only three-body force entering and Z^is also required at N 2 LO. The Ho and H2 three-body parameters are 
fixed from the 2 S 1 scattering length 133 = (0.65 + 0.04) fm and the triton binding energy B3 = 8.48 MeV, respectively. 
The Lb and Le are obtained by fixing at leading non-vanishing order by the thermal neutron-proton radiative capture 
and deuteron photodisintegration cross sections, respectively 022I1 . 

In order to compute a solution, one introduces a cutoff A which is nonphysical and thus not to be confused with 
the breakdown scale A^ of the EFT. We show cutoff variation of cross section between A = 150 MeV and A = 500 
MeV as a function of the center of mass momentum is shown in fig. 3. A numerical investigation of the Faddeev 
equation confirms these findings. One can see smooth slope and does however change significantly from order to 
order, because the dominant correction is (y/A)" and for these energy range at every order. These variations are 
nearly independent of variation of momentum. We confirm that, the cutoff variation decreases steadily as we increase 
the order of the calculation and it is of the order of (k/A)' 1 , (y/A)", where n is the order of the calculation and A = 150 
MeV is the smallest cutoff used. The other errors due to increasing momentum also appear in our calculation in low 
energy but these errors decrease when the order of calculation is increased up to N 2 LO. 




Figure 3: Curve of the cutoff variation of cross section up to N 2 LO is shown between A = 150 MeV and A = 500 MeV. The short dashed, long 
dashed and solid line correspond to LO, NLO and N 2 LO, respectively. 



The breakdown of pionless EFT in calculation is estimated to be corresponds to photon energies of m 2 n I M * 20 
MeV. The photodisintegration reaction y 3 H — » nd results in the low photon energy range(0-20 MeV) are shown in 
table 1 . It is well known that this reaction will be occurred predominantly via electric and magnetic, specially electric 
dipole transition. Our obtained results for electric dipole transition (E\) show higher contribution in comparison with 
magnetic dipole transition (Mi) at these energies. The (Ei) contribution is calculated for LO and the next leading 
order is beyond N 2 LO. The (Mi) contributions are calculated for LO, NLO and N 2 LO. The calculated cross section 
shows sharply rising from threshold to maximum about 0.88 mb at ~13 MeV and decreasing only slightly to about 
0.81 mb at ~19 MeV, in corresponding to experimental results 0]]. 

In order to provide information on the dependence of the cross section calculation on the choice of models, forces 
and currents the results for the total two-body photodisintegration cross section of 3 H have been displayed at two 
energies in table 2. At 12 MeV, one can see a good agreement between the results of AV18 and UrblX with three- 
nucleon forces and pionless EFT result up to N 2 LO. It is also found, the two-body photodisintegration of the triton 
cross section is dominated by E\ transition. 

Comparison of the different models, show that the results of pionless EFT up to N 2 LO are in good agreement with 
the other modern nucleon-nucleon potentials with three-nucleon forces. Table 2 illustrates also that the result with 
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Ey 

I 


cr(Mi) 


o-(Ei) 


cr(Mi + Ei) 


o-(Mi) 


a(Mi + Ei) 


o-(Mi) 


cr(Mi + Ei) 


Exp. 


(MeV) 


(LO) 


(LO) 


(LO) 


(NLO) 


(NLO) 


(N 2 LO) 


(N 2 LO) 


(mb) 


11.54 


0.039 


0.782 


0.821 


0.054 


0.833 


0.063 


0.845 


0.84 ± 0.05 


11.92 


0.037 


0.800 


0.837 


0.050 


0.855 


0.060 


0.860 


0.95 ± 0.05 


12.55 


0.034 


0.820 


0.854 


0.045 


0.870 


0.055 


0.875 


0.96 ± 0.06 


13.04 


0.030 


0.830 


0.860 


0.041 


0.876 


0.049 


0.880 


0.86 + 0.05 


13.92 


0.027 


0.831 


0.858 


0.038 


0.874 


0.047 


0.878 


0.79 ± 0.06 


14.92 


0.026 


0.827 


0.853 


0.036 


0.868 


0.045 


0.872 


0.71 ±0.06 


15.92 


0.023 


0.820 


0.843 


0.033 


0.858 


0.042 


0.862 


0.61 ±0.05 


16.80 


0.020 


0.810 


0.830 


0.030 


0.845 


0.038 


0.848 


0.85 ± 0.06 


17.66 


0.018 


0.796 


0.814 


0.027 


0.828 


0.035 


0.835 


0.74 ± 0.07 


18.67 


0.015 


0.786 


0.801 


0.024 


0.815 


0.032 


0.818 


0.66 ± 0.07 


19.16 


0.014 


0.778 


0.792 


0.022 


0.805 


0.030 


0.808 


0.66 ± 0.07 


19.72 


0.013 


0.770 


0.783 


0.021 


0.796 


0.029 


0.799 


0.71 ±0.11 



Table 1: The cross section for y 3 H —> nd in millibarns computed in pionless EFT up to N 2 LO. The pionless EFT cross section is comprised of the 
Ei amplitude computed to LO and the Ml amplitude computed to LO, NLO and N 2 LO. The last column shows the experimental date from QJ. 



Type of Formulation 


E y = 12 [MeV] 


E r = 40 [MeV] 


AV18-Siegert 


1.056 


0.168 


AV18-MEC 


0.949 


0.155 


CD Bonn2000-Siegert 


0.980 


0.169 


AV 1 8 +UrbanaIX-Siegert 


0.882 


0.180 


AV18+UrbanaIX-MEC 


0.915 


0.169 


CDBonn2000+TM' -Siegert 


0.889 


0.176 


EFT(N 2 LO) 


0.874 


0.240 



Table 2: The total cross section results(in mb) for two-body photodisintegration of 3 H of different models-dependent based on Faddeev ap- 
proach 1 27]. The last row shows model-independent EFT result up to N 2 LO. 
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explicit MEC differs rather substantially from that with the one-body current and implicit MEC via Siegert's theorem 
as well as EFT results. However, precise data would be required in order to reach a better conclusion. 
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Figure 4: 3 H photo-dissociation cross section <x (mb) for center of mass energy E (MeV). Dot, dash and solid curve are LO, NLO, and N 2 LO Mj 
transitions, and long-dash is LO E\ transition. 



The leading contributions from each transitions are shown in fig. 4. One can see that the M\ transition is small in 
comparisons to E\ transitions and as a results the E% transitions can be ignored for calculation in the energy range of 
our interest. A comparison between the low energy cross section computed with pionless EFT and the experimental 
data as well as the potential model calculation, can be seen in fig. 5. A detailed and very illuminating discussion of the 
experiments that contributed to these data points can be found in ref. j lj. One finds a rather good agreement between 
pionless EFT theory and experiment at low energy. Since the experimental situation is not settled, a good definite 
comparison between theory and experiment cannot be made. It is due to the fact that the two-body break up data are 
rather scattered and there is also limited experimental information on the three-body break up. 

We also compared the obtained total 3 H cross section with the other modern nucleon-nucleon AV18/UrbanaIX 
1 2711 . in fig. 5. One finds a rather good agreement between EFT results up to N 2 LO and the experimental data, as 
well as, potential models that the three-nucleon forces improve their consistency with the experiments. On the other 
hand, due to the insufficient precision of the experimental data a more conclusive comparison between theory and 
experiment cannot be made. As another result, we leave the consistency of our results for energies higher than ~ 20 
MeV. This is in agreement with the breakdown energy scale of pionless EFT (~ 20 MeV). Effects of others multipoles, 
retardation, tensor correlations, and subnuclear currents should have more and more influence at higher energies. 

Our result using EFT is model independent and universal, any model with the same input must, within the accuracy 
of our calculation, lead to the same result. Due to the large error bars in experimental data, it is not possible to draw 
further conclusions. 



4. Summery and conclusion 

We calculated the cross section of photo-dissociation of y 3 H — > nd. We applied EFT to find numerical results 
for this process. The results are compared with, AVI 8 together with the Urbana IX three-nucleon forces or CD Bonn 
with TM three-body forces, high precision other nuclear force models. Our results are in good agreement with the 
other calculations and experimental data. Theoretical predictions show the maxima sharply rising from threshold to 
maximum about 0.88 mb at ~ 13 MeV and decreasing only slightly to about 0.81 mb at ~19 MeV, in corresponding 
to experimental results. 
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Figure 5: The cross section for y^H — > nd near threshold, as a function of the incident photon energy in MeV. The solid curve corresponds to the 
cross section computed in pionless EFT up to N 2 LO. The dotes with error bar correspond to the experimental date from |l|] and triangle points 
show the Faddeev calculation results for cross sections of the 3 H photodisintegration 12711 . The short dashed, long dashed and solid line correspond 
to LO, NLO and N"LO effective field theory calculation, respectively. 



The results are quantitatively supported by available data and converge order by order in low energy expansion. 
The results are also cutoff independent when the order of calculation is increased up to N 2 LO. High quality experi- 
mental data would be very helpful more strongly to challenge theory at low energies. 
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